Coulomb Explosion and Thermal Spikes 
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A fast ion penetrating a solid creates a track of excitations. This can produce displacements 
seen as an etched track, a process initially used to detect energetic particles but now used to alter 
materials. From the seminal papers by Fleischer et al. [1] to the present [2], 'Coulomb explosion' and 
thermal spike models are treated as conflicting models for describing ion track effects. Here molecular 
dynamics simulations of electronic-sputtering, a surface manifestation of ion track formation, show 
that 'Coulomb explosion' produces a 'heat' spike so that these are early and late aspects of the same 
process. Therefore, differences in scaling are due to the use of incomplete spike models. 

PACS numbers:6f.80.Az,79.20.-m,34.50.Fa,79.20.Ap 
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Since Thompson and Rutherford there has been inter- 
est in the track of excitations produced by an energetic 
ion penetrating a solid. The fast ion excites the electron 
cloud producing excitons and electron-hole pairs in an in- 
sulator H . The electrons expelled from the track create 
additional excitations, then cool to the lattice, screening 
and finally rccombining with the holes. They eventually 
return to the ground state by radiative or non-radiative 
processes. Prior to this decay an excited region with a 
net repulsive energy persists. Fleischer et al. |j] (here- 
after FPW) proposed that this repulsion could produce 
displacements altering the material along the particle's 
track. When accounting only for the holes the process 
has been called 'Coulomb explosion', suggested to pro- 
duce amorphized tracks Q), cratering, and sputtering 
jlfl . Tracks have also been observed in liquids || , semi- 
conductors § and even metals when the excitation 
density in the track is sufficiently large. 

In describing thresholds for track registration in insula- 
tors, FPW compared their 'Coulomb explosion' model to 
the often used thermal spike model for defect production 
Such a comparison is not only of historical inter- 
est but persists in the recent literature [Q. However, one 
model describes a mechanism for energy-input (electro- 
static potential energy into kinetic energy), whereas the 
other describes the transport of energy out of a 'heated' 
cylindrical region. Here we present a molecular dynamics 
(MD) description of the response of a solid to a repul- 
sive track. We show that at high excitation densities the 
primary effect is the production of a cylindrical 'heat' 
spike. Therefore, distinctions made between 'Coulomb 
explosion' and spike models can be artifacts of the ap- 
proximations used. 

MD simulations are carried out for a model solid with 
a surface because of our interest in electronically-induced 
sputtering, the surface manifestation of track formation 
in insulators |Tc[| . As in track registration, both spike 
models [[ll) and 'Coulomb explosion' models PJT^] have 
been employed to parametrize laboratory data for elec- 
tronic sputtering. In spike models the atoms in the ion 



track are assumed to have a radial temperature pro- 
file determined by the energy deposited per unit path 
length, (dE/dx). When the energy of an atom exceeds 
a barrier, defect production in the bulk or evaporation 
at the surface can occur. From the time dependent sur- 
face temperature, T sur f, the spike sputtering yield, Y$, 
the number of atoms ejected per ion track is calculated 
integrating the surface flux over time and area. When 
a radial diffusion equation is used to determine T sur f, 
then Y s oc (dE/d x)1 ff for a fixed track radius at high 
{dE / dx) e f f [Ojrjj , where (dE/dx) e f / is the fraction of 
dE/dx going into energetic non-radiative processes [^). 
This has been used to describe sputtering produced by a 
track 'heated' by individual repulsive decay events [11|, 
secondary electrons || and a repulsive track ||Jl5],16|. 
However, when conditions leading to the (dE j dx) e ^ de- 
pendence should apply, the model fails because the en- 
ergy transport is not diffusive [ p"7| . A melt and a pressure 
pulse control the energy transport leading to [l4| ], 



Y s tt0.18(r cyl /U)(dE/dx) eff 



(0.1) 



where U is the solid's cohesive energy and the effect depth 
scales with r cy i, the initial mean radius of the spike. Here 
we show that a repulsively-induced 'heat' spike is pro- 
duced in an ion track over a broad range of excitation 
densities and the transport is not diffusive. Therefore, a 
correct spike model must be used. 

We simulate the evolution of a track of repulsive en- 
ergy that might be produced by a fast incident ion. Since 
sputtering [fl4| and crater formation in both amorphous 
and crystalline, atomic or molecular solids can be scaled 
over a broad range of material properties, we use a model 
solid made of atoms interacting via Lennard- Jones (LJ) 
potentials. We use parameters for condensed gas solids, 
but the results scale with the energy and length parame- 
ters. Also, we showed earlier that scaling was maintained 
using more complex potentials [Q. Interactions occur 
between all atoms within a cut-off radius r cut ~ 2.54Z, 
where I is the mean atom spacing. The (001) layer spac- 
ing in the fee lattice is l s « 0.8L MD simulations focusing 
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on damage have been made for repulsive energy between 
a few charges in LiF |Q , a distribution of holes in Si [jl9| , 
and using an assumed velocity distribution |Q . Here we 
follow the conversion of a track of repulsive energy into 
atomic motion in a solid. 

On electronically exciting a condensed gas solid, the 
track of excitations can produce large sputtering yields 
at high dE/dx, a process relevant to icy outer solar 
system bodies @. Ejection is presumed to be due to 
the net repulsive energy, but a quantitative description 
is lacking. At high ionization densities the net repul- 
sion can persist if the holes have low mobility and are 
not sufficiently screened during the time displacements 
could occur. Even when holes are neutralized, the ex- 
cited atoms have overlapping charged clouds which act 
repulsively at high dE/dx [0,pl|. Here we describe the 
net repulsion between 'excited' atoms in a track using 
V = (e 2 /r) exp (~r/a), where a is an average screen- 
ing constant Q . We do not distinguish between closely 
spaced, partially screened holes or overlapping excited 
species, since in both cases the electrons screen the inter- 
actions between 'excited' neighbors | |l5| , |l6[ |. The interac- 
tions of 'excited' atoms with unexcited neighbors was left 
unchanged; including a weak polarization had a small ef- 
fect. For energy conservation, a larger cut-off, r™"' = 7a, 
is used for interactions between 'excited' species. 
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FIG. 1. (a). The sputtering yield from a repul- 
sive track vs. [l a (dj/dx)]. The solid lines are 
Yfu — 14. In [1.3 (a/r CO ui)] [Is dJ/dx] 2 . [Results scale with 
LJ parameters. Those used here are e = 10.3 meV and 
a = 3.405 A giving a cohesive energy characteristic of a 
number of condensed-gas solids, U = 0.08 eV and l s = 2.66 
A.] (b). Yield scaled to the fit Y$u in Fig. la vs. r, the 
neutralization/recombination time. These are roughly fit by 
exp[— a(TD/r) x ] with a and x varying slowly with [l s dj/dx]; 
curve shown is for a = 1.13i s , [hdJ/dx] = 1 with a — 0.46 
and x = 0.56. 



The material is 'excited' by instantaneously changing 
the potentials between atoms in the track. The resulting 
velocities and positions of all particles are then followed 
p7[ . This is done for a number of excitation densities 
and screening constants. In a second set of calculations 
the excitations were quenched (recombination) . This was 
done statistically so the average number of excitations in 
the track decayed exponentially, exp [— t/r] where r is the 
quenching (neutralization) time. The number of excita- 
tions per unit path length, dJ/dx, is related to dE/dx 
for fast ions |IJ and the mean radius of the distribu- 
tions, r cou i, depends on the speed of the incident ion. 
For given track parameters [dJ/dx, r cou i (here ~ l s )] ex- 
cited species are chosen randomly. On excitation, atoms 
are sputtered if they cross a plane 2r™"' above the 'top', 
here the (001) surface. The number of atoms ejected 
in each run is called the yield, which is related to the 
track of damage in FPW but is more easily measured. 
Sample depth was chosen to be at least twice r£°"'. Dif- 
ferent boundary conditions did not change the average 
yields significantly. The sample size (3a; 10 4 — 3xl0 5 ) 
and simulation times (15 — 80 ps) were adjusted to a 
and dJ/dx. Extending times by tens of picoseconds or 
doubling the thickness did not change the average yield. 
The distribution in the yield is broad, especially for small 
l s (dJ/dx) where sputtering occurs when two excitations 
are produced close together [|H]. Results, averaged over 
~ 10—200 excitation distributions, are given in Fig. la as 
a function of charge density, l s (dJ/dx), for 3 values of the 
screening constant a. The yield is seen to be quadratic in 
dJ/dx for each a and increases non-linearly with a; lines 
Y = 14. In [1.3 (a/ w)] [h {dJ/dx)f give a good fit [|. 
The yields as a function of neutralization time, r, for a 
fixed (dJ/dx) and a are given in Fig. lb scaled by the 
fit in Fig la. Because neutralization is treated stochas- 
tically, the size of the yield is affected even for relatively 
large r and only becomes independent of the neutraliza- 
tion rate for t > 10t£>, with tjj the Debye period for the 
lattice (~ 0.5 ps) |jl7). However, the yield is still seen to 
be nearly quadratic in (dJ/dx) for the full range of r with 
an additional weak dependence on excitation density. 

In MD simulations details of the energy transfer and 
sputtering can be extracted. At very low excitation den- 
sities (Y < 1) the probability of neighbors being excited 
near the surface leads to ejection. At the highest excita- 
tion densities (Y > 5), the repulsive energy is transferred 
to neighbors in times ~ 0.2td producing a cylindrically 
heated region. Comparing the evolution of the result- 
ing radial temperature profiles with our earlier simula- 
tions of 'heat' spikes, the energy transport is found to 
be similar. That is, we find for the repulsively produced 
heat spike transport is not diffusive Wj\ , a melt and a 
pressure pulse are formed and these control the energy 
transport. Seiberling et al. JT2I suggested that a repulsive 
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track could produce a heat spike but used the standard 
diffusive model to obtaining a yield dependence very dif- 
ferent from that calculated here. 
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FIG. 2. Energy distribution of the ejecta for [l s dj/dx] — 2, 
a = 1.13Z S and r cou i — l s , for t — td (a) and r = 20td 
(b). The dotted line shows the energy spectrum from a cylin- 
drical spike, as in ref. (uj, with (dE/dx) eff « 45U/L and 
r cy i ~ 2.6 Is- For the "Coulomb" spike produced in (b), 
(dE/dx) e ff ~ 38U/l s and r cy i m 2.6 l s , while the spike in (b) 
has a 30% lower (dE/dx) e ff. Peaks are prompt ejecta, with 
position determined by the potential energy between neigh- 
bors and the surface binding energy. 

From the energy distributions for the ejecta at high 
excitation density two regions can be seen in Fig. 2. At 
large ejecta energies, E > 10[/, the spectra has peaks 
due to prompt ejection from the upper surface layers of 
the initial track. This accounts for ~ 20% of the ejecta 
at large dJ/dx but dominates at very small dJ/dx. On 
the other hand, the principal component of the ejecta 
in Fig. 2 exhibits an energy distribution like that found 
in our studies of ejection from a narrow cylindrical heat 
spike shown as the dashed line [jlTj. That is, there is a 
broad, quasi-thermal distribution at low ejecta energies, 
E < U which gives way to a non-Maxwellian, ~ E~ 2 
dependence at E > U. The latter dependence is charac- 
teristic of low energy cascades in a solid [Q but differs 
from thermal spike model predictions. Surprisingly, al- 
though the ejection process changes in going from the 
lowest to the highest dJ/dx, there is no dramatic change 
in the dependence of the yield on dJ/dx in Fig. 1. This 
is due to the nearly linear dependence of spike yield on 
(dE/dx) e ff discussed below. 

From the above, the primary effect of the repulsive 
decay at the high excitation densities is the production 
of a 'heat' spike. This is also shown quantitatively us- 
ing the MD calculation to determine the repulsive en- 
ergy driving atomic motion in times short compared to 
desorption times. At ~ 0.2t£> the energy density de- 
posited via repulsion is (dE / dx) rep w 0.l5e 2 (dJ/dx) 2 for 
a = I.13l s and r cou i w l s . This energy is localized in a 
cylindrical region of radius r rep w 2.6 l s for our values of 
r CO ui- Using (dE/dx) e ff = (dE/dx) rep and r cyl = r rep 
as initial conditions in the expression for the yield from 



a cylindrical 'heat' spike in Eq. |0.l| gives a yield that has 
the same dependence on dJ/dx and a size within 10% 
of those in Fig. la. Therefore, if a repulsive region is 
sustained, a 'heat spike' is formed at high dJ/dx which 
determines the subsequent energy transport, sputtering 
and displacements. Therefore, 'Coulomb explosion' and 
spikes are the early and late aspects of the same pro- 
cess, so that contradictions for effects occurring later 
than ~ 0.2td are due to incorrect descriptions of the 
'heat' spike. Below we examine under what conditions a 
'Coulomb' track is sustained long enough to produce a 
'heat' spike. 
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FIG. 3. The change in the total energy in the track of inter- 
acting holes and electrons vs. time for 'excited' electrons with 
[Is dJ/dx] = 1, and a dielectric constant e — 1.6e - A fraction 
of the energy will be converted to atomic motion. No scat- 
tering or inelastic energy loss (solid), no inelastic energy loss 
but scattering X 3 — I (dashed), for \i n — I andAe; n = 0.03 
eV (dotted), for A^ n = I and Aei„ = 0.1 eV (dash-dotted). 

Since neutralization quenches the energy of a spike, we 
simulated the neutralization and screening of holes by 
free electrons. That is. we make a separate simulation 
in which we track the cooling of the 'excited' electrons 
in the field of the track of positive charges. The lattice 
is like that in the MD simulations and we choose appro- 
priate mean-free paths for elastic, A s , or inelastic, A-j n , 
scattering, both assumed to be isotropic |22| ]. A time de- 
pendent a or an average a and t can be determined from 
such a simulation. Instead, we give in Fig. 3 the instan- 
taneous total potential energy density in the track vs. 
time. This is the energy available for repulsive heating of 
the lattice. The interactions between the positive charges 
in the track and the electrons are Coulombic beyond an 
atomic radius, a . Inside a the potentials gradually be- 
come flat giving a binding energy of ~ 8.6 eV. We ex- 
amined a number of potential forms inside a Q and initial 
electron energies and the trends are the same. Electrons 
are 'excited' by receiving, on the average, 15 eV, roughly 
the energy dissipated to the lattice (W-value minus ion- 
ization energy |22|] ). This does not treat the very fast 
electrons which would further slow neutralization. 

For only elastic scattering, the total energy in the track 
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in Fig. 3 is sustained with a very large effective a. Ex- 
amining the radial dependence of the electron density, a 
fraction of the electrons cool rapidly by electron-electron 
collisions partially neutralizing the track. This is di- 
electronic recombination, also called Auger recombina- 
tion. The remainder of the electron cloud, now 'hotter' 
but more fully screened, expands to a larger average ra- 
dius. Varying A s affects the electron density distribu- 
tion, whereas including an inelastic energy loss, Aej n , 
decreases the potential energy in the track. This also 
heats the lattice, but here we are interested in the repul- 
sive energy. For a Aei n and Ai„ roughly corresponding 
to electrons in a polymer or in ice (dotted curve) [ ^2f , it 
is seen in Fig. 3 that most of the initial track energy re- 
mains during the time it takes to produce a 'heat spike' 
repulsively, ~ 0.2td. This would correspond to an in- 
termediate t and relatively large a in Figs, lb and 2a. 
Therefore, for a reasonable inelastic energy loss, screen- 
ing is not sufficient to quench the repulsive production 
of a 'heat spike'. Increasing Ae;„ (dot-dashed curve) or 
decreasing \ n leads to more rapid but still incomplete 
quenching over the relevant time. Quantities for specific 
measurements need to be used and the role of excited 
states in the track must be included. 

In this paper we close the circle on one aspect of an 
old [Q but topical problem, the effect on a solid of a 
track of excitations produced by a fast penetrating ion. 
A MD simulation of repulsive explosion in a track of ex- 
citations has been carried out and the effect of neutral- 
ization has been examined. This was done for narrow 
screened Coulomb spikes in a condensed-gas solid to con- 
strain the simulation size. However, we showed earlier 
that the yields scale with U and [dE/dx) e tf. Whereas 
at low excitation densities ejection can occur if neighbor- 
ing 'excited' species are formed near the surface HQ, at 
the higher excitation densities the repulsive energy in the 
track produces a 'heat' spike. Therefore, 'Coulomb ex- 
plosion' and spike models are the early and late aspects 
of the repulsive decay process. The spike formed by the 
repulsive heating of the lattice can lead to displacements 
and ejection (electronic sputtering), examined here. Be- 
cause the yield is roughly linear in (dE / dx) e f / at high 
excitation densities and (dE/dx) e ff is proportional to 
(dJ/dx) 2 , the yield from a repulsive track is seen (Fig. la) 
to vary smoothly, having the same quadratic dependence 
on dJ/dx in going from low to high excitation densities al- 
though the character of the ejection process changes. At 
high excitation densities the yield can now be estimated 
by substituting the repulsive energy deposition into the 
new, useful expression for spike sputtering in Eq. p.l| . 
At higher excitation densities than those studied here, 
the more distributed 'heat' spike produced by the elec- 
tron cooling to the lattice might also be important ||. 
Earlier applications of a 'heat' spike for describing track 
formation or electronic sputtering foundered on the use 
of a spike model in which the energy transport was de- 



scribed incorrectly. That is, the transport is not diffusive 
]l4] , p7t and the local kinetic energy distribution is not 
Maxwellian (Fig. 2). Calculations should now focus on 
accurately describing the interactions in an ion track for 
times < 0.2 to including electrons, holes and excitons, 
as the latter also contribute to the net repulsive energy. 
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